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RECENT  DEVELOPMENTS  IN  THE 
THEORY  OF  ATMOSPHERIC  ROTORS 

by  James  D.  Doyle  and  Dale  R.  Durran 


The  Sierra  Nevada  Range  is 
one  of  the  most  prominent 
and  steepest  mountain  barri¬ 
ers  in  the  United  States  and,  not 
surprisingly,  is  a  well-known  lo¬ 
cation  for  a  multitude  of  topo¬ 
graphically  forced  atmospheric 
phenomena.  As  the  prevailing 
westerly  winds  pass  over  the  Si¬ 
erra  Nevada,  gravity  waves  are 
frequently  generated.  Occasion¬ 
ally  these  mountain  waves  result 
in  spectacular  topographically 
forced  phenomena  such  as 
trapped  lee  waves,  downslope 
windstorms,  rotors,  and  atten¬ 
dant  wave  and  rotor  cloud  struc¬ 
tures,  as  shown  in  the  examples  in 
Fig.  1.  In  situations  such  as  those 
shown  in  Fig.  1,  severe  downslope 
winds  near  the  surface,  some¬ 
times  in  excess  of  50  m  s-1,  decel¬ 
erate  rapidly  in  the  lee  and  give 
way  to  a  return  flow  back  toward 
the  mountain  crest  that  is  the 
lower  branch  of  an  intense  hori¬ 
zontal  circulation,  as  illustrated  in 
the  schematic  in  Fig.  2.  These 
horizontal  vortices,  known  as  ro¬ 
tors,  are  common  to  the  steep 
eastern  slopes  of  the  Sierra  Nevada,  and  particularly 
over  the  Owens  Valley,  where  they  are  notorious  for 
their  intensity  (e.g.,  Whelan  2000).  Rotors  have  also 


Fig.  I.  Photographs  over  the  Owens  Valley  in  the  lee  of  the  Sierra  Ne¬ 
vada  Range,  taken  during  the  Sierra  Wave  Project,  illustrating  (a)  and  (b) 
common  rotor  char¬ 
acteristics  such  as  ro- 
tor  and  lenticular 
clouds  and  (a)  blow¬ 
ing  dust.  The  flow  in 
both  photographs  is 
from  right  to  left  as 
viewed  from  the 
north. 
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Fig.  2.  Schematic  streamlines  illustrating  a  rotor  circu¬ 
lation  and  attendant  cloud  features  (adapted  from 
Ludlam  and  Scorer  1 957).  Regions  of  clear-air  turbu¬ 
lence  associated  with  the  rotor  circulation  are  denoted 
by  the  red  symbols. 


been  observed  in  a  number  of  other  mountainous 
regions,  including  the  Rockies  (Lester  and  Fingerhut 
1974)  and  various  locations  in  Europe  (e.g.,  Queney 
et  al.  1960).  Rotors  can  be  severe  aeronautical  hazards 
and  have  been  cited  as  contributing  to  numerous  air¬ 
craft  encounters  with  severe  turbulence  and  accidents, 
including  occurrences  involving  modern  commercial 
and  military  aircraft.  Rotor  circulations  may  also  be 
important  for  the  lifting  and  transport  of  aerosols  and 
contaminants,  as  apparent  in  Fig.  la. 

Mountain  waves  and  rotors  were  the  subject  of  one 
of  the  first  modern  U.S.  multiagency  field  programs 
in  meteorology,  the  Sierra  Wave  Project  (Holmboe 
and  Klieforth  1957;  Grubisifand  Lewis  2004,  personal 
communication)  that  took  place  in  the  early  1950s 
(during  which  the  photographs  in  Fig.  1  were  taken). 
Although  research  on  mountain  waves  and  related 
phenomena  has  been  very  active  in  the  50  years  since 
the  Sierra  Wave  Project,  little  attention  has  been  de¬ 
voted  to  rotors.  Furthermore,  there  have  been  almost 
no  thorough  observations  of  real-world  rotors  because 
they  are  dangerous  to  sample  using  in  situ  aircraft 
measurements,  and  their  small  spatial  scale  prevents 
them  from  being  resolved  by  conventional  observing 
networks.  Therefore,  it  is  not  surprising  that  moun¬ 
tain-induced  rotors  still  remain  poorly  understood 
and  difficult  to  forecast.  Only  recently  have  compu¬ 
tational  resources  been  sufficient  to  enable  what  are 
thought  to  be  the  first  realistic  numerical  simulations 
of  rotors.  Here  we  present  results  from  a  high-reso¬ 
lution  nonhydrostatic  mesoscale  model,  the  Coupled 


Ocean-Atmosphere  Mesoscale  Prediction  System 
(COAMPS™;  Hodur  1997),  which  provides  a  fresh 
perspective  on  the  structure  and  characteristics  of 
mountain-induced  rotors. 

Numerical  simulations  have  suggested  that  a  key 
aspect  of  rotor  development  involves  the  synergistic 
interaction  between  the  lee  wave  and  boundary  layer 
effects  as  discussed  recently  in  Doyle  and  Durran 
(2002).  As  an  example,  consider  the  two-dimensional 
numerical  simulation  of  flow  over  bell-shaped  topog¬ 
raphy  shown  in  Fig.  3,  which  is  initialized  based  on  a 
blend  of  the  potential  temperature  and  cross-moun¬ 
tain  wind  components  observed  in  the  Grand  Junc¬ 
tion,  Colorado;  Denver,  Colorado;  and  Lander,  Wyo¬ 
ming,  soundings  at  1200  UTC  3  March  1991.  On  this 
day,  strong  winds  were  observed  on  the  lee  side  of  the 
Colorado  Front  Range  and  rotor  clouds  were  reported 
at  the  U.S.  Air  Force  Academy  in  Colorado  Springs, 
Colorado,  between  1355  and  1555  UTC.  The  horizon¬ 
tal  resolution  of  these  simulations  is  100  m  with  95 
vertical  levels  spanning  the  interval  0  <  z  <  11.6  km. 
The  details  of  the  model  configuration  are  described 
by  Doyle  and  Durran  (2002).  After  3  h  of  model  in¬ 
tegration,  a  series  of  trapped  lee  waves  has  developed 
in  the  streamline  field  with  trough-to-crest  ampli- 


Fig.  3.  Streamlines  and  horizontal  vorticity  (/-compo¬ 
nent)  from  a  two-dimensional  numerical  simulation 
initialized  based  on  a  sounding  representative  of  the 
conditions  in  central  Colorado  at  1 200  UTC  3  Mar  1991. 
The  cross-mountain  wind  speed  less  than  or  equal  to 
zero  is  shown  using  blue  isotachs  (every  I  m  s_l).  The 
horizontal  vorticity  greater  than  0.02  s_l  is  shaded  in 
color.  (From  Doyle  and  Durran  2002.) 
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tudes  exceeding  1200  m.  Intense  downslope  winds, 
with  a  maximum  of  25  m  s_1,  develop  along  the  lee  slope. 
Regions  of  reversed  cross-mountain  winds  associated 
with  a  series  of  rotors  are  present  near  the  surface  be¬ 
neath  the  crests  of  the  first  three  waves.  The  first  rotor 
downstream  of  the  ridgeline  is  associated  with  a  maxi¬ 
mum  reversed  cross-mountain  wind  component  of 
5.9  m  s-1  and  a  region  of  recirculating  flow  that  extends 
above  the  mountain  crest.  A  sheet  ofy-component  vor- 
ticity  ( du/dx  -  dw/dz)  originates  in  the  region  of  high 
shear  within  the  boundary  layer,  which  forms  due  to 
surface  friction  processes  along  the  lee  slope.  The  vor¬ 
tex  sheet  separates  from  the  surface,  ascends  into  the 
crest  of  the  first  lee  wave,  and  remains  aloft  as  it  is  ad- 
vected  downstream  by  the  undulating  flow  in  the  lee 
waves.  Boundary  layer  separation  is  facilitated  by  the  ad¬ 
verse  pressure  gradients  associated  with  trapped  moun¬ 
tain  lee  waves.  Vertical  and  horizontal  wind  shear  and 
turbulence  production  are  maximized  along  the  el¬ 
evated  sheet  of  horizontal  vorticity,  particularly  along 
the  leading  edge  of  the  rotor  circulation  (as  shown 
schematically  in  Fig.  2),  in  general  agreement  with 
anecdotal  evidence  from  the  Sierra  Wave  Project 
(Holmboe  and  Klieforth  1957)  and  aircraft  observa¬ 
tions  of  rotors  taken  in  the  lee  of  the  Rocky  Mountain 
Front  Range  analyzed  by  Lester  and  Fingerhut  (1974). 

Boundary-layer  processes  are  of  fundamental  im¬ 
portance  in  determining  rotor  development  and  in¬ 
tensity.  For  example,  an  otherwise  identical  simula¬ 
tion  to  that  shown  in  Fig.  3,  in  which  a  free-slip  lower 
boundary  condition  was  applied,  failed  to  produce 
rotors  or  reversed  flow.  Transient  rotors  can  be  gen¬ 
erated  in  the  presence  of  a  free-slip  lower  boundary 
using  a  highly  idealized  upstream  sounding,  but  re¬ 
alistic  rotors  with  characteristics  similar  to  those 
found  in  nature  appear  to  develop  only  in  the  pres¬ 
ence  of  surface  friction.  A  series  of  experiments  con¬ 
ducted  by  Doyle  and  Durran  (2002)  indicate  that  an 
increase  in  the  surface  roughness  beyond  a  small  value 
weakens  the  reversed  flow  and  decreases  the  depth  of 
the  rotor  circulation.  The  surface  roughness  also  has 
an  influence  on  the  downstream  location  of  the  rotor 
because  of  the  impact  of  the  lee  wave  pressure  per¬ 
turbations  on  the  boundary  layer  separation.  An 
asymmetry  in  the  solar  radiation  reaching  the  ground 
commonly  occurs  in  mountain-wave  events  because  of 
the  presence  of  foehn  clouds  upstream,  while  the  lee 
slope  remains  predominantly  cloud  free.  An  additional 
set  of  simulations  was  conducted  that  indicate  increas¬ 
ing  the  lee-side  heat  flux  deepens  the  rotor  circula¬ 
tion,  in  general  agreement  with  observations  taken 


during  the  Sierra  Wave  Project  (Kuettner  1959).  The 
turbulence  intensity  within  the  rotor  also  increases 
due  to  the  contribution  of  heating  from  the  surface. 

In  order  to  examine  the  dependence  of  the  rotor 
characteristics  on  the  lee-wave  amplitude,  Doyle  and 
Durran  (2002)  conducted  an  additional  series  of 
simulations  that  were  performed  with  an  atmospheric 
structure  in  which  the  upstream  static  stability  is  con¬ 
stant  throughout  each  of  two  horizontal  layers.  These 
simulations  used  varying  mountain  heights  and  inter¬ 
face  depths  to  examine  a  portion  of  the  parameter 
space  relevant  to  basic  rotor  dynamics.  The  results 
indicate  that  the  magnitude  of  the  reversed  flow  in  the 
primary  rotor  for  a  simulation  with  surface  friction 
is  highly  correlated  with  the  strength  of  the  adverse 
pressure  gradient  in  the  primary  lee  wave  in  an  oth¬ 
erwise  identical  simulation  without  surface  friction 
(i.e.,  with  a  free-slip  lower  boundary  condition). 
These  simulations  also  show  that,  even  with  surface 
friction,  no  rotors  form  unless  the  adverse  pressure 
gradient  in  the  corresponding  free-slip  simulation 
exceeds  a  threshold  value,  further  demonstrating  syn¬ 
ergistic  coupling  between  boundary  layer  induced 
vorticity  and  lee-wave-induced  adverse  pressure  gra¬ 
dients  in  the  formation  of  low-level  rotors. 


Fig.  4.  GOES  visible  satellite  image  valid  at  0000  UTC 
29  October  2000. 
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Fig.  5.  SkewT/log-P  thermodynamic  sounding  for 
Vandenberg,  CA,  for  0000  UTC  29  October  2000.  The 
horizontal  wind  speed  (m  s_l  )  and  direction  are  shown 
to  the  right.  One  full  wind  barb  corresponds  to  5  m  s_l 
and  a  flag  denotes  25  m  s_l. 


signature  of  a  single  stationary  wave  to  the  south.  A 
sounding  from  Vandenberg,  California,  located  ap¬ 
proximately  150  km  upstream  of  the  Sierra  Nevada, 
at  0000  UTC  29  October  2000  is  shown  in  Fig.  5.  The 
sounding  and  the  larger-scale  COAMPS  simulation 
suggest  that  strong  cross-mountain  flow  (10-15  m  s_1 
at  700  hPa  and  20-25  m  s_1  at  500  hPa)  and  an  inver¬ 
sion  near  the  elevation  of  the  Sierra  crest  existed  up¬ 
stream  of  the  mountains,  characteristics  often  asso¬ 
ciated  with  severe  downslope  winds  along  the  lee 
slopes  (e.g.,  Durran  1990)  and  rotor  circulations. 

The  Sierra  Nevada  is  a  continuous  and  narrow 
mountain  range  with  a  mean  width  of  approximately 
100  km  and  a  length  of  650  km.  The  geographical  re¬ 
gion  of  interest  for  the  numerical  simulation  is  fo¬ 
cused  on  the  Owens  Valley,  California,  which  is  ad¬ 
jacent  to  the  topography  comprised  of  the  main  crest 
of  the  high  Sierra  to  the  west,  with  peaks  over  4400  m 
elevation,  that  abruptly  descends  to  the  valley  floor  at 


High-resolution  simulations  of 
rotors  indicate  the  presence  of 
small-scale  vortices  or  subrotors 
embedded  within  the  parent  ro¬ 
tor.  These  subrotors,  which  in 
some  respects  are  analogous  to 
multiple  suction  vortices  within  a 
tornado,  may  pose  the  greatest 
hazard  to  aviation.  Vortex 
stretching  and  tilting  are  likely  to 
play  an  important  role  in  the  en¬ 
hancement  of  vorticity  within 
such  subrotors,  but  these  pro¬ 
cesses  can  play  no  role  in  chang¬ 
ing  the  vorticity  in  a  two-dimen¬ 
sional  flow,  such  as  that  shown  in 
Fig.  3. 

To  further  examine  the  inter¬ 
nal  structure  of  rotors  found  in 
nature,  a  three-dimensional  high- 
resolution  real  data  simulation 
was  carried  out  for  a  case  of 
mountain  waves  on  28  and  29 
October  2000  in  the  lee  of  the  Si¬ 
erra  Nevada  range.  The  GOES 
visible  satellite  image  for  0000 
UTC  29  October  2000,  shown  in 
Fig.  4,  indicates  numerous 
trapped  mountain  waves  in  the 
lee  of  the  northern  Sierra  and  a 
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Fig.  6.  Topography  (m)  for  the  fourth  grid  mesh  (Ax  =  I  km)  is  shown 
by  the  color  scale  and  contours  every  250  m.  The  location  of  the  fifth 
grid  mesh  (Ax  =  333  m)  is  represented  by  the  white  rectangle.  The 
dashed  latitude  and  longitude  lines  are  shown  every  0.25°. 
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Fig.  7.  Streamlines  and  horizontal  vorticity  (y- component)  from  a  three-dimensional  numerical  simulation  at 
the  1 2-h  time  valid  at  0000  UTC  29  Oct  2000  for  (a)  along  axis  AB  (see  Fig.  6)  and  (b)  approximately  3  km  south 
of  axis  AB.  The  horizontal  vorticity  greater  than  0.005  s_l  is  shaded  in  color.  The  cross-mountain  wind  speed  less 
than  or  equal  to  zero  is  shown  using  blue  isotachs  (every  I  m  s_l). 


a  mean  elevation  of  1200  m,  as  shown  in  Fig.  6.  The 
eastern  portion  of  the  valley  is  bounded  by  the  White 
and  Inyo  Mountains.  Near  Independence,  California, 
the  width  of  the  Owens  Valley  is  approximately  15  km. 

The  COAMPS  model  simulation  is  initialized  at 
1200  UTC  28  October  2000  using  all  available  obser¬ 
vations  blended  with  a  first  guess  based  on  the  Navy 
Operational  Global  Atmospheric  Prediction  System, 
which  is  also  used  for  lateral  boundary  conditions. 
In  this  application,  the  model  makes  use  of  five  nested 
grid  meshes  centered  on  the  Owens  Valley,  with  the 
finest  mesh  using  a  333-m  grid  increment  (location 
shown  in  Fig.  6).  A  west-east  oriented  vertical  cross 
section  of  streamlines  and  horizontal  vorticity  along 
line  AB  (Fig.  6)  at  the  12-h  time  valid  at  0000  UTC 
29  October  2000  is  shown  in  Fig.  7a.  Similar  to  the 
idealized  two-dimensional  simulation  shown  in 
Fig.  3,  a  sheet  of  large  horizontal  vorticity  (y-com- 
ponent)  is  significantly  enhanced  in  the  boundary 
layer  along  the  lee  slopes  of  the  Sierra  range  and  lifted 
upward  into  the  lee  wave  crest  to  near  7  km.  Two 
small-scale  circulations  are  present  within  the  rotor 
near  the  surface,  perhaps  indicative  of  the  role  of 
three-dimensionality  and  the  complex  forcing  asso¬ 
ciated  with  the  upstream  and  downstream  topogra¬ 
phy.  The  strong  horizontal  and  vertical  shears,  par¬ 
ticularly  at  the  upstream  edge  of  the  rotor,  generate 
intense  turbulence  that  is  advected  into  the  upstream 
edge  of  the  lee  wave. 


An  indication  of  the  along-ridge  variation  in  the 
finescale  structure  within  the  rotor  may  be  obtained 
by  examining  a  second  vertical  cross  section  shown 
in  Fig.  7b,  which  is  located  approximately  3  km  to  the 
south  of  that  shown  in  Fig.  7a.  The  flow  at  4  km  in 
Fig.  7b  is  strongly  reversed,  with  maximum  speeds  of 
5  m  s-1  that  are  more  intense  than  the  extremum  in 
the  reversed  flow  near  the  surface.  The  substantial 
variation  in  rotor  substructure  between  the  cross  sec¬ 
tions  in  Figs.  7a, b  is  consistent  with  the  possibility  that 
vortex  stretching  and  tilting  may  play  a  role  in  am¬ 
plifying  any  subrotors  that  do  develop.  The  substruc¬ 
tures  shown  in  Figs.  7a, b  are  also  highly  transient,  sug¬ 
gesting  that  a  given  area  within  the  rotor  can  be 
relatively  benign  at  one  moment  and  extremely  dan¬ 
gerous  a  few  minutes  later.  Such  localized  regions  of 
turbulence  and  shear  within  the  parent  rotor  may  be 
one  of  the  key  terrain-induced  phenomena  that  con¬ 
tribute  to  an  enhanced  rate  of  aviation  accidents  in 
mountainous  regions  (e.g.,  Carney  et.  al  1995). 

Guidance  for  forecasting  rotors  is  virtually  nonex¬ 
istent  despite  their  important  and  direct  link  to  avia¬ 
tion  safety  issues.  The  general  characteristics  of  rotors 
have  been  known  since  the  completion  of  the  Sierra 
Wave  Project  in  the  1950s,  yet  surprisingly  little  is 
known  about  several  of  their  key  aspects,  including 
their  internal  structure,  strength,  climatology,  and  mi¬ 
crophysical  processes.  Recent  numerical  simulations 
suggest  that  the  meteorological  conditions  favorable 
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for  the  formation  of  strong  rotors  are  similar  to  those 
that  favor  the  development  of  trapped  lee  waves  and 
downslope  winds,  namely  1)  a  significant  cross¬ 
mountain  wind  component  of  at  least  10  m  s-1;  2)  an 
elevated  inversion  near  crest  level;  and  3)  a  significant 
increase  in  the  cross-mountain  wind  speed  above  crest 
level.  Nevertheless,  more  extensive  theoretical  and  ob¬ 
servational  efforts  are  required  to  verify  and  refine 
these  general  guidelines.  Forecasting  of  mountain-wave 
activity  and  associated  turbulence  in  support  of  aviation 
operations  are  commonly  performed  using  empirical 
and  statistical  techniques,  which  could  be  extended  to 
include  rotors  if  more  sophisticated  guidelines  prove  to 
be  reliable. 

In  the  five  decades  following  the  Sierra  Wave 
Project,  the  characteristics  and  dynamics  of  rotors 
have  remained  an  enigma.  New  observations  of  ro¬ 
tors,  such  as  from  the  forthcoming  Terrain-Induced 
Rotor  Experiment  (T-REX;  Grubisifand  Kuettner 
2003),  will  be  needed  to  gain  insight  into  the  charac¬ 
teristics  of  rotors,  including  their  internal  structure, 
and  to  evaluate  the  results  from  high-resolution  and 
eddy- resolving  numerical  models  to  determine  the 
predictability  of  rotors.  The  multiphase  T-REX  field 
campaign  and  high-resolution  numerical  model 
simulations  should  ultimately  provide  the  means  to 
substantially  improve  our  understanding  and  ability 
to  forecast  mountain-wave  induced  rotors. 
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